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Abstract. We exploit the correspondence between the three-dimensional Lorentzian 
Einstein-Weyl geometries of the hyper-CR type, and the Veronese webs to show that 
the former structures are locally given in terms of solutions to the dispersionless Hirota 
equation. We also demonstrate how to construct hyper-CR Einstein-Weyl structures by 
Kodaira deformations of the flat twistor space TCP 1 , and how show to recover the pen- 
cil of Poisson structures in five dimensions illustrating the method by an example of the 
Veronese web on the Heisenberg group. 



1. Introduction 

The notion of three-dimensional Veronese webs appearing in the study of finite dimen- 
sional bi-Hamiltonian system is based on existence of one-parameter families of foliations of 
an open set of K 3 by surfaces [151 Ell EZ] • The same structure underlies the Einstein-Weyl 
geometry in 2+1 dimensions, where the surfaces are required to be totally geodesies with 
respect to some torsion-free connection compatible with a conformal structure [61 [16] . The 
connection between the Veronese webs and Einstein-Weyl geometry has not so far been 
made, and one purpose of this note is to show that the three-dimensional Veronese webs 
correspond to a subclass of Einstein-Weyl structures which arise as symmetry reductions 
of (para) hyper-Hermitian structures in signature (2, 2). This class of Einstein-Weyl struc- 
tures is called hyper-CR as it admits a hyperboloid of (para) CR structures |11| . In the 
next Section we shall make use of this correspondence to show that all hyper-CR Einstein- 
Weyl structures locally arise from solutions of the dispersionless Hirota equation (Theorem 
12. ip . In Section [3] we shall elucidate the procedure of recovering a hyper-CR Einstein-Weyl 
structure (and thus also a Veronese web) from the corresponding twistor space. In Section 
H] (Theorem I4.1|) we give an explicit construction of the bi-Hamiltonian structure in five 
dimensions from the hyper-CR Einstein-Weyl structures. In this construction the twistor 
function corresponds to the Casimir of the Poisson pencil. Finally we use the example of 
a Veronese web on the three-dimensional Heisenberg group to illustrate the procedure of 
recovering a solution of the dispersionless Hirota equation from a given three-parameter 
family of twistor curves. 



2. Einstein-Weyl structures from dispersionless Hirota equation 

Consider the dispersionless integrable equation of Hirota type [311 [22] (see also [8] for 
a discussion of its discrete version) 

(b - a) w x w yz + a w y w zx - b w z w xy = 0, (1) 

where w = w(x, y, z) is the unknown function on an open set BcR 3 and a, b are non-zero 
constants such that a ^ b. This equation arises as the Frobenius integrability condition 
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[Lq,Li] = for the dispersionless Lax pair of vector fields 

L = d z d x + A a d z , L x = d y - —d x + A b d y . (2) 

w x w x 

This Lax pair is linear in the spectral parameter A and does not contain derivatives with 
respect to this parameter. It therefore fits into the formalism of [11], and one expects 
equation ([1]) to give rise to an Einstein- Weyl structure on B of hyper-CR type. To construct 
this structure we could find a linear combination of Lq,Li of the form V\ — AV2,V2 — 
AV3, where Vi,i = 1,2,3 are vector fields on B, then read off the contravariant conformal 
structure V2 V2 — V\ V3 and try to solve a system of differential equations for the 
Einstein- Weyl one-form. We shall instead use another procedure which will give the one- 
form directly (yet another, more straightforward method applicable to a broad class of 
dispersionless integrable PDEs has been proposed recently [14j). Extend ([2]) to a Lax pair 
of vector fields on a four-manifold M = BxR, where r is the coordinate on M 

Lo' = Lo + d T , L\ = L\. 

This Lax pair is of the form L ' = W l - \W 2 ,L X ' = W 3 - AW 4 , where Wi, W 2 , W 3 , W4 
are linearly independent vector fields on M. Therefore it defines an anti-self-dual hyper- 
hermitian conformal structure of neutral signature |10l [12], [7] given in its contravariant form 
b^O W\ W4 — W 2 W3. Let g be the corresponding covariant form. This hyper-hermitian 
structure admits a non-null Killing vector K = d/dr, and thus the Jones-Tod procedure 
|17j gives rise to an Einstein- Weyl structure on the space of orbits B of the group action 
generated by K in M. This is in general given by 

h= \K\- 2 g- \K\- 4 K®K, uj = 2\K\~ 2 * g (K A dK), (3) 

where K = g(K, •), \K\ 2 = g(K, K), and * g is the Hodge operator of g. This defines a Weyl 
structure consisting of a conformal structure [h] = {ch,c : B — > and a torsion-free 

connection D which is compatible with [h] is the sense that 

Dh = uj h. 

This compatibility condition is invariant under the transformation h — > 4> 2 h, uj — > lo + 
2<i(ln0), where is a non-zero function on B. 

Applying the Jones-Tod procedure ([3]) to the Lax pair ([2]) , rescaling the resulting three- 
metric by cj) 2 = w z /(w x w y ), and using ([1]) to simplify the resulting one-form yields the 
Weyl structure 

h = -^d X 2 + ( d y 2 + ( ^ d* (4) 

w y w z [a - by w x w z (a - b) z w x w y 

a b ab 

+2- dxdy - 2- — dxdz + 2- --3 dydz, 

[a — b) w z [a — b) w y (a — b) z w x 

to = dx — ay dz. 

w x w„ w z 



^The vector fields L'§,L' X span an RP 1 worth of surfaces - called a surfaces - through each point of 
M. The resulting conformal structure [g] is uniquely defined by demanding that these surfaces are totally 
isotropic and self-dual. The Frobenius integrability condition [Lo,^'i] = modulo L' Q and L[ is equivalent 
to vanishing of the self-dual part of the Weyl tensor of [g]. In general the distribution Z/o,£'i i s defined 
on the RP 1 bundle of real projective spinors over M, and the vector fields L' Q ,L'i contain derivatives w.r.t 
the parameter A. If no such derivatives are present, then [g] contains a metric which is hyper-Hermitian 
|10l 112] . i. e. there exist three real anti-commuting endomorphisms I,S,T : TM — > TM such that 

S 2 =T 2 = 1ST = -I 2 = 1, 

and any g € [g] is hermitian w.r.t the hyperboloid a 2 — b 2 — c 2 = 1 of complex structures J = al + bS + cT. 
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A Weyl structure is said to be Einstein-Weyl if the symmetrised Ricci tensor of D is 
proportional to some metric h € [h]. This conformally invariant condition can be formulated 
directly as a set of non-linear PDEs on the pair (h,co): 

Rij + 2 + ~\{ R+ \^ kt °k + j^ k ^kj hij = 0, 

where V, Rij, and R are respectively the Levi-Civita connection, the Ricci tensor and the 
Ricci scalar of h. A direct computation shows that the Weyl structure is Einstein- 
Weyl iff the dispersionless Hirota equation (TjQ) holds. We must have w x w y w z ^ for 
the conformal structure to be non-degenerate. For example w = x + y + z gives a flat 
Einstein-Weyl structure, with oj = and h the Minkowski metric on B = K 2 ' . 

According to Gelfand and Zakharevich [IS], a Veronese web on a three-manifold B is a 
one-parameter family of foliations of B by surfaces, such that the normal vector fields to 
these surfaces through any p £ B form a Veronese curve in F(T*B). The Veronese web 
underlying the PDE (TjQ) corresponds to a dual Veronese curve in ¥(T p B) 

A — ► V(A) = Vi- 2XV 2 + AV 3 , 
where A £ RP , and the vector fields V{ are given by 

IV ID 

V\ = ab(aw y d z - bw z d y ), V 2 = ab(w y d z - w z d y ), V 3 = bw y d z - aw z d y + (a - b) y z d x . 

One of the seminal results of [31J is that all Veronese webs in three dimensions are locally 
of this form, and thus arise from solutions to equation (TIJ. 

We observe that h(V(X), V(A)) = for all A, where h is given by ([2]), so V(A) determines 
a pencil of null vectors. The vector fields V\ — XV 2 and V 2 — AV3 (or equivalently the 
vector fields Lq and L\ given by ^)) form and orthogonal complement of V(A). For 
each A € MP 1 they span a null surface in B which is totally geodesic with respect to 
D. The existence of such totally geodesic surfaces for a triple (B, [h],D) is equivalent to 
the Einstein-Weyl condition [6|. Thus Veronese webs in three dimensions form a subset 
of Einstein-Weyl structures. The necessary and sufficient condition for an Einstein-Weyl 
structure to correspond to a Veronese web is that its Lax pair does not contain derivatives 
w.r.t. A. We have therefore established 

Theorem 2.1. There is a one to one correspondence between three-dimensional Veronese 
webs and Lorentzian Einstein-Weyl structure of hyper-CR type. All hyper-CR Einstein- 
Weyl structures are locally of the form where the function w = w(x,y,z) satisfies 
equation ([!]). 

This Theorem, together with results of [TU] implies that there is a one to one correspon- 
dence between hyper-CR Einstein-Weyl structures and third order ODEs such that certain 
two fibre preserving contact invariants vanish. 

Twistor theory. In the real analytic category one can make use of a complexified setting, 
where w is a holomorphic function on a complex domain Be C C 3 , and the Einstein-Weyl 
structure is also holomorphic. The twistor space of the PDE (TIJ is a two-dimensional 
complex manifold Z whose points correspond to totally geodesic surfaces in Bq. This gives 
rise to a double fibration picture 

B c £ B c x CP 1 A Z, (5) 

where Z arises as a factor space of Be x CP 1 by the rank-two distribution {Lq, Li} spanned 
by ([2])- The points in Bq correspond to rational curves, called twistor curves, in Z with 
self-intersection two, i.e. the normal bundle is N(l p ) = 0(2), where l p = \x o v~ l {p), for 
P G B c . 
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The special case of the hyper-CR Einstein- Weyl structures is characterised by the ex- 
istence of a holomorphic fibration of Z over a complex projective line The existence 
of this fibration is a consequence of the absence of vertical djd\ terms in the Lax pair 
([2]). The real manifold B (and thus the real solutions to ([T|)) corresponds to twistor curves 
invariant under an anti-holomorphic involution r : Z Z which fixes an equator of each 
twistor curve. 

Hierarchies. Equation (pQ) can be embedded into an infinite hierarchy of integrable PDEs, 
which arises from the distribution 

Li = di - T7—d x + A ai di, i = 0, 1, . . . 
o x w 

and is given by [Li, Lj] = 0, or 

(oj — aj) d x wdidjW + ajdiwdjd x w — aidjwdid x w = 0, (no summation). (6) 

Here di = d/dx l and cjj are distinct constants. The function w = w(x,xo,x±, . . . ) depends 
on an infinite number of independent variables, but the Cauchy data for the overdetermined 
system d5J) only depends on functions of two-variables which is the same functional degrees 
of freedom as in ([Tj). 

3. Hyper-CR equation and Kodaira deformation theory 

In [11] it was shown that all Lorentzian hyper-CR Einstein-Weyl spaces on an open set 
Bel 3 are locally of the form 

h = (dY + H x dT) 2 - A(dX - H Y dT)dT, u = H XX dY + {H x H X x + 2H XY )dT, (7) 

where H = H(X, Y, T) satisfies the the hyper-CR equation EH [13 [261 ED S] 

Hxt — Hyy + HyHxx — H x H X y = 0. (8) 

Combining this result with Theorem 12. II suggests that PDEs ([1]) and (jHJ) are equivalent, as 
there exists a local diffeomorphism of B mapping Q to ([7|). It can be shown that there is 
no point equivalence between these equations |5J, but this does not rule out the existence 
of a Backlund type transformation between them which involves the second derivatives of 
the dependent variables (a transformation of this type connects the first and the second 
Plebanski heavenly equations |12j ) or a Legendre transformation of the type recently anal- 
ysed by Bogdanov [2]. A procedure recovering a solution of the Hirota equation JT]) from 
the hyper-CR equation (JSJ) will be illustrated (at the end of Section by the example of 
the Veronese web on the Heisenberg group. 

In this section we shall simplify the procedure of recovering the solution to ([8]) and 
the corresponding Einstein-Weyl structure form the twistor data Equation ([8]) is 

equivalent to [Lo,£i] = 0, where the Lax pair ia^ 

L = d Y -X(d T + H Y dx), h = d x - \(dy + H x d x ). (9) 
The ring of twistor functions spanning the kernel of this Lax pair is (ip, A), where 

^ = T + XY + X 2 X + \ 3 H + A 4 V> 4 + • • • = Xi ( 10 ) 

i=0 

and the functions tpi depend on (X, Y, T) but not on A. The recursion relations connecting 
ipi+i to tpi arise from equating coefficients of A in Lqi^ = Liip = to zero. The consistency 
conditions for these relations imply that any of the functions ipi satisfies 

(d x d T -d Y + B Y d\ - H x d x dy)A = 0. 



This differs from the Lax pair of [TT] by replacing A by 1/A. 
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Note that this is different than a linearisation of ([5]). 

In the real analytic setup, where the twistor methods can be applied, the conformal 
structure [h] arises from defining a null vector at p € Be to be a section of a normal bundle 
N(l p ) which vanishes at one point on the twistor line l p = CP 1 in the twistor space Z to 
second order. To define the connection D, define a direction at p € Be to be the one- 
dimensional space of sections on N(l p ) which vanish at two points. The one-dimensional 
family of curves vanishing at these two points corresponds to a geodesies in Be- In the 
limiting case when these two points coincide this geodesic is null w.r.t [h] in agreement with 
the definition of the Weyl structure [16]. The procedure of recovering both D and [h] from 
a family of curves in Z is explicit, but rather involved [23]. This can be vastly simplified 
by using the local form ([7]): given a three-parameter family of sections of Z — > CP 1 

A— ►(A,^(A)) 

choose a point on the base of the fibration Z — > CP 1 , say A = 0, and parametrise this 
family by coordinates (X, Y, T) where 

dip 1 d 2 ip 

T = ip\ x=0 , Y = —\ x=0 , X = --^\ x=0 . (11) 

Expanding the twistor function ip in a Laurent series in A gives ()10p . and we can read off 
H = H(X,Y,T) from the coefficient of A 3 . The Einstein- Weyl structure is now given by 
, and we did not have to follow the procedure of [23] . 

For example, expanding the Nil twistor function |11] (which depends on a constant 
parameter e, and in the limit e — > reduces to a section A — > T + A Y + X 2 X of the 
undeformed twistor space TCP 1 ) 

^ = T + XY In (1 - XeX) (12) 

in a power series in A and comparing with (|1U|) gives H = eX 2 /2. Using ([7D we find the 
corresponding one-parameter family of deformations of the flat Einstein- Weyl structure is 
given by Lorentzian Nil geometry 

h = (dY + eXdT) 2 - AdXdT, uj = e(dY + eX dT) . (13) 

We shall come back to this example at the end of Section HI 

Deformation theory. Let [tp : ttq : iri] be homogeneous holomorphic coordinates on the 
twistor space Z = TCP 1 corresponding to the flat Einstein-Weyl structure with H = 0. 
Here [ttq : tt\] are homogeneous coordinates on the base CP 1 and A = Cover Z 

by two open sets U and U with 7ri ^ on U and ttq ^ on U. Following the steps of 
Penrose's non-linear graviton construction |29] we construct twistor spaces corresponding 
to non-trivial Einstein-Weyl spaces by deforming the patching relation between U and hi. 
The Kodaira theorem [18] guarantees that the deformations preserve the three-parameter 
family of curves and the deformed twistor space still gives rise to a three-dimensional 
manifold (Be, [h],D) with Einstein-Weyl structure. 

The hyper-CR deformations preserve the fibration of Z over CP 1 . Thus we deform the 
patching relations such that ttq/tti is preserved on the overlap. The infinitesimal deforma- 
tions of this kind are generated by elements of -ff 1 (-Z, 0), where denotes a sheaf of germs 
of holomorphic vector fields. Let a vector field 

„ d ( d d \ 

J = J a7 + 9\ 7T °a h7Fl ^ — ) 

Olp V OTTq OTTi J 
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defined on an overlap of U and U define a class in H 1 ^, 0). Here the functions / and g 
are holomorphic and homogeneous of degree two and zero respectively on the overlap. The 
finite deformation ip = ip(ip, A, e) is given by integrating 

dip d-iro diri 

de de de 

It is possible that the assumption about real-analyticity can be dropped by employing 
the techniques of [M] or some modification of the dispersionless 8 approach [3j, but that 
remains to be seen. 



4. Bi-Hamiltonian structure 

Recall that a Poisson structure on a manifold U is a bivector P € A 2 (TU) such that the 
associated Poisson bracket {/,<?} := P-i (df A dg) satisfies the Jacobi identity 

{f,{g,h}} + {h,{f,g}} + {g,{h,f}} = 

for all smooth functions f,g, h on U. If we regard U C M. n as an open set, and adapt local 
coordinates such that 

P = y P ^A_ A JL 

a, 13=1 

then the Poisson bracket takes the form 

yQ /3 df dg 



{f,g}= £ P 



dx a dxP ' 

a,/3=l 

If n is odd, the Poisson structure necessarily admits at least one Casimir, i.e. a non-constant 
function C such that {C, /} = for all functions / on U. 

A bi-Hamiltonian structure P x on U consists of two Poisson structures Pq and P\ such 
that the pencil of Poisson brackets 

{f,g}x ■= {f,g}o + Hf,g}i 

satisfies the Jacobi identity for all values of the parameter A (see e.g. [20J [23 [1] ) . We shall 
assume that P^ has exactly one Casimir function C A for each A. Thus, regarding the P x € 
r(A 2 (T[/)) as a map from T*U to TU, we find that this map has a one-dimensional kernel. 
Therefore, for each A G MP 1 , the manifold U is foliated by four-dimensional symplectic 
leaves of P A , i. e. there exists a rank four integrable distribution T>\ C TU such that 
T>\ = ker (dC x ). This distribution is annihilated by a A-dependent one-form 

e(A) := {P x AP A )Jft, 

where f2 is some fixed volume form on U, and J denotes its contraction with a section of 
A 4 (TU). Therefore e(A) is quadratic in A, and we can write 

e(A) = e 3 + Ae 2 + AV (14) 

for some one-forms e\ The Pfaff theorem implies that dC x G span(e(A)). 

The intersection of all distributions T>\ as A € MP 1 varies is an integrable rank-two 
distribution which we shall call T>. Thus 

V = ker e 1 n ker e 2 n ker e 3 = H ker(e(A)) C TU 

A 

which is integrable, as an intersection of tangent bundles to the symplectic leaves of the 
Casimir C A . Therefore V C V x C TU. We shall assume that all these distributions 
have constant ranks. The two dimensional distribution T>z '■= T>\/T> is defined on the three 
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dimensional quotient B := U/T>. This distribution is integrable, and it defines the structure 
of a Veronese web on B. This web is a projection of symplectic leaves of P x to B. 

This is essentially the construction of [IS]. In [TS] an inverse construction of a bi- 
Hamiltonian structure from a Veronese web is also given. This construction does not 
appear to be explicit, but it has been shown [3D] that a bi-Hamiltonian structure on U can 
be recovered form a Veronese web on B. 

We shall now give a simple algorithm for recovering the bi-Hamiltonian structure from 
any hyperCR Einstein- Weyl structure which, by Theorem 12.11 is equivalent to a Veronese 
web. In our procedure the distribution T>z is identified with the twistor distribution ([2]) 
given by the span of the Lax pair Lq and L\ which (in the complexified setting) underlies 
the double fibration picture (J5]). The one forms e % in (|14[) define the conformal structure 

h = e 2 ® e 2 - 2(e 1 ®e 3 + e 3 ® e 1 ) 

and the Probenius condition e A de = is equivalent to the hyper-CR Einstein-Weyl 
equations [11] . 

Before formulating the next Theorem recall that the conformal signature of signature 
(2, 1) on a three-dimensional manifold B is equivalent to the existence of an isomorphism 

TB^SQS, (15) 

where S is a rank two real symplectic vector bundle over B. In concrete terms (|15p identifies 
vectors with two by two symmetric matrices (or equivalently, with homogeneous quadratic 
polynomials in two variables). The conformal structure is then defined by declaring a vector 
to be null iff the corresponding matrix has rank one (or equivalently, if the corresponding 
homogeneous polynomial has a repeated root). The isomorphism (|15|) is a particular exam- 
ple of a GL(2, R) structure: the group GL(2, H.) acts on the fibres of S, and induces Lorentz 
rotations and conformal rescalings on vectors in B as 0[(2, M) = so(2, 1) © R. 
Let (po,Pi) be local coordinates on the fibres of the dual vector bundle §*. 

Theorem 4.1. Let {B, [h],D) be a hyper-CR Einstein-Weyl structure and let (Lq,Li) be 
the commuting twistor distribution (Lax pair) spanning a pencil of null and totally geodesic 
surfaces through each point p G B. Then 

P X = L A — + L 1 A— . (16) 
opo dpi 

is a bi-Hamiltonian structure on the five-dimensional manifold §*, where S* is a real rank 
two bundle over B such that H5\) is the canonical isomorphism given by the conformal 
structure [h] on B. 

Proof. Consider three linearly independent vector fields (Vi,V2,V3) on B such that 
the conformal structure in the Einstein-Weyl structure is represented by the contravariant 
form V2 V2 — V\ © V3. Let (\ C B be a two-dimensional totally geodesic surface through 
p € B corresponding to A G MP 1 and let (po,pi) be coordinates on the fibres of the 
cotangent bundle T*(\. The tangent space T p Q\ is spanned by the twistor distribution 
Lq = V\ — XV2,Li = V2 — AV3 (any twistor distribution underlying a hyper-CR Einstein- 
Weyl space can be put in this form, by taking a linear combination of the spanning vector 
with A-independent coefficients) and so the canonical Poisson structure on T*£a is 

d r d 

^oA — +Li A— . 
opo opi 

We now extend this to a bi-Hamiltonian structure on BxM 2 . To do it we need to canonically 
identify the tangent spaces to different null totally geodesic surfaces in B and show that 
the union of cotangent bundles Uaprp 1 ^*Ca is isomorphic to S*. To establish the latter, 
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observe that under the isomorphism (|15|) vectors tangent to £a correspond to quadratic 
polynomials with one of the roots equal to A: represent the basis (Vx, V2, V3) of T p B as a 
symmetric matrix of vector fields Vab with A, B = 0, 1 and Vbo = Ui, Vqi = Vio = V2, V\\ = 



B 



are homogeneous 



V3, and the twistor distribution is La = k^Vab where 7r 

coordinates on MP 1 such that ir 1 /ir° = —A. Thus any vector tangent to Cx is W = \i a La 
for some [i A = \pp,fj}\. The two by two matrices corresponding to such vectors are W 



7T, with ir A fixed and fi A varying. Equivalently vectors tangent to (\ correspond to 
homogeneous polynomials with common root Dividing out by the factor containing 

this root yields linear homogeneous polynomials, i.e. an isomorphism T(\ = S|a . To obtain 
the canonical identification of all tangent spaces T p Ca as A varies we simply evaluate the 
twistor distribution at different values of A: Let £1 and (2 be totally geodesic null surfaces 
corresponding to Ai and A2 respectively. The isomorphism 0X2 : PCi ~~ > ^2 is defined on 
the basis vectors by 

<^12(-^0 1 Ai) = Lq\\ 2 , <Al2(^lUi) = Li\\ 2 . 

Therefore the fibre coordinates (po,pi) can be used unambiguously for all A and the Poisson 
structure on T*(\ extends to a Poisson pencil (fT6|) on Uagirp 1 ^*Ica = 

For any Poisson structure given by a bivector P the Jacobi identity reduces to 



5=1 



. n. 



This set of conditions holds, with n = 5, for P x = P = Pq — APi regardless of a particular 
form of the Lax pair (Lq,Li), as long as [Lo,^i] = because in this case the Frobenius 
theorem can be used to put the Poisson bi-vector in the canonical form. 

If the vector fields of the Lax pair are linear in the parameter A - which is the case for 
Lax pairs underlying hyperCR Einstein- Weyl structures - then the bivector (|16p defines a 
pencil of compatible Poisson structures. 

□ 



In the construction of the Poisson pencil 
we called ip in formula (|10p . and 



TBI) the Casimir C A is a twistor function which 



T^z = span(Lo, Li), T> = span(<9/<9po> 9/ dpi). 
Using the explicit form of the twistor distribution ([2]) yields the following 

Corollary 4.2. Let w = w(x,y,z) be a function on B with continuous second derivatives. 
The Poisson brackets 



Pn 



( 



w z /w x 







-1 

-1 





1 






-Wy/W X \ 

1 




/ 



Pi 



/o 





V 










define a compatible Poisson pencil if and only if the function w 
dispersionless Hirota equation ([I]). 



w(x,y,z) 



\ 
b 



°/ 

satisfies the 



^This agrees with the Gelfand-Zakharevich construction [TS], where the existence of the isomorphism 
T*B = S* S* is postulated. The authors consider 

T* B ® §* = (§* 0S*0 §*) (§* ® A 2 (§*)). 

The 5-manifold supporting the Poisson pencil arises as the second factor in this decomposition, i.e. the total 
space of §* — > B as the line bundle A 2 (S*) is trivialised. 
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To obtain a bi-Hamiltonian system, the Hamiltonian function needs to be selected as one 
of the coefficients in the expansion of the Casimir C A in the powers of A. It follows from 
the general theory that all coefficients of this expansion are in involution with respect to 
both Poisson structures Pq and P\. Thus the coordinate functions on B are first integrals 
of the Hamiltonian flow, and the integral curves of this flow are vertical to B. They are of 
the form Pk{t) = Uk(x,y, z)t + pk(x,y,z), where k = 0,1, and Uk,Pk are known functions 
on B whose form depends on the choice of the Hamiltonian. This simple form of solution 
curves reflects the fact that the Gelfand-Zakharevich construction of Veronese webs puts 
the system in the action-angle coordinates. 

Example. The Heisenberg group. We have already shown that the Heisenberg group 
carries an Einstein- Weyl structure of Lorentzian signature given by (|13p . The twistor 
distribution ([9]) defining the Veronese web is 

L = dY-Xd T , Li = dx - A(<9y + eXdx)- (17) 

The metric h defining the conformal structure is left-invariant with the Killing vectors 
generated by the right invariant vector fields 

R Y = dy, R T = dr, R x = d x - eTdy 

forming the Heisenberg Lie algebra 

[R x ,R T ] = eR Y , [R x ,R Y ]=0, [R t ,R y } = 0. 

The Einstein- Weyl structure is only invariant under the abelian subalgebra spanned by Ry 
and Rt as the Lie derivative of the twistor distribution along Rx is non-zero. The Poisson 
pencil is also not invariant under the Heisenberg action, which is best seen by writing it 
in terms of the left-invariant vector fields (which commute with the right-invariant vector 
fields) as 

P X = (L Y A 8 po + Lx A d Pl ) — \{{L T + eXL Y ) A d po + (L Y + eXL x ) A d pi ) 

and C Rx P x = eX P° / 0. We note that there is one Poisson bracket in the pencil (corre- 
sponding to A = 0) which is left-invariant. 

We finish the discussion of this example by presenting the solution to the dispersionless 
Hirota equation (HD corresponding to the Heisenberg Veronese web (|17|) . The procedure we 
shall use elucidates two different parametrisations of twistor curves leading to equations (JSJ) 
and (P) respectively. In case of the hyper-CR equation we choose a point, say A = 0, on 
the base of the fibration Z — > CP 1 (here we again work in the complexified setting), and 
parametrise the curve l p given by A — > (A, ip{\)) corresponding to p G B by the position of 
its intersection with the fibre over A = and the first two jet coordinates as in The 
third jet then defines the dependent variable H in equation ([5]). This can be seen from 
formula (fT2j) . 

In case of the dispersionless Hirota equation we instead choose four distinct points 
Ai,A2,A3 and A4 on the base. The intersections of the curve with the three fibres over 
Ai,A2,A3 give, up to a reparametrisation, the coordinates (x,y,z). The intersection with 
the fibre over A4 then defines the dependent variable w. We can use the Mobius transfor- 
mation to map (Ai, A2, A3) to (0, 00, 1). Thus we expect the coordinate T in the hyper-CR 
equation to map directly to coordinate x in the dispersionless Hirota equation. This is the 
twistorial interpretation of the original procedure of Zakharevich where the two-dimensional 
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distribution T>z denning the Veronese web is 

V z {\) = ker ((Ai - A 4 )(A - A 2 )(A - X 3 )w x dx 

+ (A 2 - A 4 )(A - Ai)(A - \?)wydy 

+ (A 3 -A 4 )(A-A 1 )(A-A 2 Kdz). 

Comparing this with T>z(X) = span(Lo, L\) given by the twistor distribution (JTTJ) and using 
the Mobius freedom described above yields T>z(0) = ker (dx),T>z(l) = ker (dz), T>z(oo) = 
ker (dy), or equivalently 

dijj\x=o £ span(<ix), dip\x=i € span(cte), dip\\ =00 € span(dy), 

where ip is the twistor function (|12p . These formulae define the coordinates up to reparametri- 
sations (x,y,z) — > (x(x),y(y), z(z)), and we choose 

X = T, y = Xe~ eY , z = (1 - e X)e~ e( - T+Y \ 

The solution w(x, y, z) of the dispersionless Hirota equation (pQ) can now be read-off from 
T^z{Xa) = ker (dw) which yields w = (1 — e\4 L X)e~ e ^ i lT + Y ) ( or after transforming to the 
new variables and reabsorbing some constants by rescaling (x,y,z), 

w = ye ax + ze bx , 

where A 4 = 1 — b/a. 
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